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Introduction 

Unlike CO, which is observed extensively in the inter- 
stellar medium of our own and other galaxies, its homo- 
logue SiO is observed principally in outflows associated 
with regions of star formation. This striking difference in 
behaviour is a consequence of the lower elemental abun- 
dance and the more complete depletion of silicon from the 
gas phase during grain formation. Both carbon and sili- 
con form refractory materials - graphite and silicates - 
of which the cores of interstellar grains are believed to 
be composed; but the much higher elemental abundance 
of carbon, nc/nn — 3.55 x 10~ 4 , compared with silicon, 
nsi/n-n = 3.37 x 1CT 5 (Anders & Grevesse 1989), leads to 
some of the carbon remaining in the gas phase. 

The SiO molecule was first detected in the Galactic cen- 
tre (Sgr B2) by Wilson et al. (1971) and subsequently in Ori 
A by Dickinson (1972). More recent observations of SiO in 
jets (see, for example, Bachiller et al. 1991, Martin-Pintado 
et al. 1992, Codella et al. 1999, Nisini et al. 2007) imply 
that, in these objects, at least some of the silicon has been 
restored to the gas phase; this can be achieved through 
sputtering of the grain material, probably in shock waves, 
which are features of the outflows. It has been known since 



the study by Draine et al. (1983) that grain ice-mantles 
can be eroded in C-type shock waves, owing to impact of 
neutral particles on charged grains at the ion-neutral drift 
speed, which is the speed of ambipolar diffusion. More re- 
cent work (Flower & Pineau des Forets 1995, Flower et al. 
1996, Field et al. 1997, May et al. 2000) has shown that 
this process might result also in the partial erosion of the 
refractory grain cores. The simulations undertaken by May 
et al. of the sputtering of various silicates (forsterite, fay- 
alite and olivine) by neutral atoms showed that C-shock 
speeds in excess of 30 km s _1 are necessary to erode a 
significant fraction (more than a few per cent) of these ma- 
terials. On the other hand, it has since been recognized 
(Le Bourlot et al. 2002, Ciolek et al. 2004) that the speeds 
at which C-type shock waves can propagate are limited, 
both by the collisional dissociation of molecular hydrogen, 
which leads to a thermal runaway and the formation of 
a J-type shock wave, and by the ion magnetosonic speed, 
whose value is constrained by the inertia of the charged 
grains. Consequently, the efficacy of the erosion of silicates 
in C-type shock waves is restricted by the maximum C- 
shock speed, which depends on the preshock density and 
transverse magnetic field strength, and on the fraction of 
charged grains (Flower & Pineau des Forets 2003). 

In a prev ious study of SiO pro duction in the interstel- 
lar medium, Schilke et al. (1997) (henceforth Sch97) con- 
sidered in detail the erosion of silicon from grain cores and 
from their mantles by C-type shock waves and the resulting 
SiO emission spectrum; this study remains the only one of 
its kind that has been published to date. In the interven- 
ing decade, there has been progress in our understanding 
of the sputtering process (May et al. 2000), the gas-phase 
chemistry of silicon (Le Picard et al. 2001), and the grain 
dynamics (Flower & Pineau des Forets 2003), and so it 
seems timely to revisit this topic. Recent observations of 
SiO in jets (Nisini et al. 2007) extend to higher rotational 
levels than previously and provide a further motivation for 
such an study. Accordingly, we have reconsidered the chem- 
istry of silicon in steady-state, planar MHD shock waves, 
with a view to providing a grid of models which may serve 
in the analysis of current and future observations of rota- 



tional transitions of SiO in outflows. These models yield 
additional results relating to rovibrational transitions of 
H2 and rotational lines of CO, as well as forbidden lines 
of atoms and atomic ions, including [Fc II] (cf. Giannini et 
al. 2006). 

Grain-grain collisions and the sputtering of silicon 
in oblique C-type shock waves, in which the preshock 
magnetic field direction is inclined to the plane of the 
shock front, have been considered by Caselli et al. (1997). 
However, such simulations have yet to incorporate the gas- 
phase chemistry and the radiative cooling processes, in or- 
der to enable quantitative analyses of the spectral line ob- 
servations to be made. 



1. The model 

We summarize below those developments, in both the MHD 
code and the data used and produced by the code, which are 
relevant to the modelling of the intensities and the profiles 
of the rotational lines of SiO. We take as our baseline the 
study by Sch97. The reader is referred to the more recent 
papers cited below for details of the modifications. 

1.1. The dynamics of charged grains 

The main revision of the code itself concerns the treat- 
ment of the charge and of the dynamical effects of the 
grains. The charge distribution of both the grains and 
the "very small grains" (VSG), represented by poly- 
cyclic aromatic hydrocarbons in our model, are calcu- 
lated, assuming that collisions with gas-phase particles 
(electrons, ions and neutrals) dete rmine the charge; see 
Flower & Pineau des Forets (2003) As was mentioned in 
the Introduction, allowance for the mass density of the 
(mainly negatively) charged grains can reduce significantly 
the magnetosonic speed in the charged fluid 
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where T± are the temperatures, ji± are the mean masses 
and p± are the mass densities of the positively and nega- 
tively charged fluids; B is the magnetic field strength in 
the preshock gas. The magnetosonic speed is the maxi- 
mum speed at which a C-type shock can propagate in the 
medium. Furthermore, the momentum transfer between the 
charged grains and the neutral fluid affects the ion-neutral 
drift speed and has consequences for the degree of sputter- 
ing of the grains within a C-type shock wave. 

1.2. Radiative cooling by H2 

The thermal balance of the medium, particularly the ra- 
diative cooling by H2, is treated more exactly in the cur- 
rent model. Rovibrational excitation of H2, principally 
by H, H2 itself, and He, followed by radiative decay, is 
the princi pal mechanism of cooling the shock-heated gas. 
Following |Le Bourlot et al. (2002) the populations of the 
rovibrational levels of H2 are calculated in parallel with 
the hydrodynamical and chemical rate equations, yield- 
ing the most accurate determination of the rate of cooling 
by H2 that is achievable within the framework of a time- 
independent (steady-state) model of the shock structure. 
The rate coefficients for the collisional excitation by H of 
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rovibrational transitions of H2 are from the recent work of 
Wrathmall et al. (2007) Collisional dissociation and ion- 
ization of H2, as well as ionization of H, are taken into 
account. Collisional dissociation of H2 is a particularly im- 
portant process, as the removal of H2 can lead to a thermal 
runaway. The associated increase in the kinetic tempera- 
ture, T n , of the neutral fluid, and hence of the adiabatic 
sound speed 

2 5fc B T n 

where [A n is the mean mass of the neutral fluid, can give rise 
to a sonic point in the flow and hence to a J-type discon- 
tinuity. This phenomenon imposes an additional constraint 
on the maximum speed of a C-type shock wave in a molec- 
ular medium. 

In the compressed gas of the postshock region where 
most of SiO forms, cooling by 12 CO, 13 CO, and H 2 starts 
to dominate that by H2 . In order to predict accurately the 
emitted SiO emission spectrum, an LVG treatment of the 
cooling by th ese species has been introduc ed, following the 
procedures of Neufeld & Kaufman (1993) A "thermal gra- 
dient" c s /z', where z' = z + 1.0 x 10 cm and z is the 
independent integration variable, is added quadratically to 
the macroscopic velocity gradient, in order to simulate pho- 
ton escape through thermal line broadening. 



1.3. The sputtering of grains 

The sputtering probabilities computed by May et al. (2000) 
for olivine (MgFeSiC^) are used to determine the rate 
of erosion of Si from (charged) silicate grains by neutral 
particles. We include also the sputtering of carbonaceous 
(amorphous carbo n) grains, using the sputtering yields of 
Field et al. (1997) Impacts (at the ion-neutral drift speed) 
of abundant heavy neutral species are the most effective 
in eroding the grain cores. For collisions with CO, for ex- 
ample, we adopted the sputtering probabilities calculated 
for impacts of Si atoms, which have the same mass as CO 
and hence the same impact energy. As May et al. (2000) 
have shown, similar yields of Si are obtained from the three 
types of silicate: fayalite, Fe2Si04; forsterite, Mg 2 Si04j and 
olivine, MgFeSi0 4 . 

The grain mantles are eroded first, through impacts 
with the most abundant neutral species, H, H2 and He, 
at ion-neutral drift speeds which are below the thresh- 
old for erosion of the cores (Draine et al. 1983, Flower & 
Pineau des Forets 1994). The initial composition of the 
gas is calculated in chemical equilibrium, whilst that of 
the grain mantles, and the elemental depletion into the 
grain cores, is based on observations (cf. Flower & Pineau 
des Forets 2003, table 1). We incorporated a representa- 
tive polycyclic aromatic hydrocarbon (PAH), with a frac- 
tional abundance npAH/"-H = 10 -6 , an upper limit which 
is consistent with estimates of the fraction of elemental 
carbon likely to be present in the form of very small 
grains (Li & Draine 2001; Weingartner & Draine 2001). 
Th e state of charge of this species was computed follow- 
ing Flower fc Pineau des Forets (2003)| who showed that 
increasing the fractional abundance of the PAH results in 
higher values of the magnetosonic speed in the preshock 
gas, thereby enabling C-type shock waves to propagate at 
higher speeds. In their turn, the higher speed shocks erode 



the silicate grains more effectively. However, the adopted 
PAH abundance has little effect on the internal structure 
of the shock wave, as the grains become negatively charged 
in the shock wave and dominate grain-neutral momentum 
transfer. 

The total number density of the grains was computed 

— from the total mass of the refractory cores, which is 
0.78 x 10~ 2 times the mass of the gas for the adopted 
composition of the cores; 

— assuming a bulk mass density of 3 g cm~ 3 ; and 

— taking a size distribution of the grain cores 
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following Mathis et al. (1977)| where n g (a g ) is the num- 



ber density of grains with radius a g , assumed to have 
upper and lower limits of 0.3 /jm and 0.01 /j,m, respec- 
tively. 

The thickness of the grain mantles was determined from 
their molecular composition (Flower & Pineau des Forets 
2003, table 2), assuming a mean number of 5 x 10 4 molecu- 
lar binding sites per layer of the mantle and a thickness 
of 2 x 10~ 4 /im for each layer, independent of the size 
of the grain core; there are no Si-containing species in 
the mantles. This procedure yields an initial mantle thick- 
ness of 0.015 /xm on the grain cores, whose mean radius is 
a g = 0.02 /j,m. However, the erosion of the grain mantles 
occurs sufficiently rapidly, as the ion and neutral flow veloc- 
ities begin to decouple in the shock wave, that the existence 
of thick ice-mantles on the grains in the preshock gas has 
little effect on the shock dynamics (see fig. 6 of Guillet et 
al. 2007). 

1.4. Chemistry 

The chemical reaction network comprises over 900 re- 
actions connecting the abundances of more than 100 
species, in both the gas and the solid phases. The 
complete list of species and reactions is available from 



|http: / / massey.dur.ac.uk/ drf/ outflows-test /species-chemistryjshock/ 
In the context of the present study, we emphasize the 
gas-phase chemistry of silicon and the formation of SiO, 
in particular. The total rate of cosmic ray ionization of 
hydrogen, (, was taken to be 5 x 10~ 17 s _1 . 

Two reactions are important for the formation of 
SiO from Si, which is released into the gas phase by the 
sputtering of the grains, namely 



Si + 2 — ► SiO + O 
Si + OH — > SiO + H 
for which the following rate coefficients (cm 3 s _1 ) 

fci = 2.7 x 10- 10 exp(-lll/T) 
k 2 = 1.0 x 10- 10 exp(-lll/T) 



(1) 
(2) 

(3) 
(4) 



were adopted by Sch97 from the compilation of 
Langer fc Glassgold (1990)] The exponential factor in 

that 



equation Q3J) derives from the argument (Graff 
the reactions proceed only with the fine-structure states 
of Si (3p 2 3 Pj) with J > 0, of which J = 1, which lies 
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111 K above the J = ground state, is the more signifi- 
cantly populated at low temperatures. More recently, the 
rate coefficient for reaction ([I]) has been measured at low 
temperatures (15 < T < 300 K) by |Le Picard et al. (2001) 
and found to be given by 



fci = 1.72 x 10- 10 (T/300)-°- 53 exp(-17/T). (5) 

We have adopted ([5|) and the same expression for &2- 
Evidently, the differences between the present and previ- 
ous values of these rate coefficients are most significant for 
temperatures T <; 100 K, i.e. in the cooling flow of the 
shocked gas. 

The abundance of SiO is limited by its conversion to 
SiC>2 in the reaction with OH 



SiO + OH 



SiO, 



H 



(G) 



whose rate coefficient remains subject to considerable un- 
certainty. We adopt the same expression as Sch97, viz. 



k e = 1.0 x 10" n (T/300) 
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in units of cm 3 



However, we note 



(7) 
that 



Zachariah & Tsang (1995) calculated a barrier of 433 K to 
reaction ([6J), and a rate coefficient 

k 6 = 2.5 x 10" 12 (T/300) - 78 exp(-613/T); 



see the discussion of |Le Picard et al. (2001)[ At T = 300 K, 
the latter rate coefficient is 30 times smaller than the for- 
mer. In the ambient (preshock) and the postshock gas, 
where T w 10 K, the existence of an activation energy of 
several hundred kelvin would prevent the oxidation of SiO 
in reaction ([6]) from occurring. The rate coefficient for re- 
action © in the UMIST data base (Le Teuff et al. 2000) 
is 

fc 6 = 2.0 x 10~ 12 

in cm 3 s _1 , which is 50 times smaller than ([7|) at T = 10 K. 
Fortunately, the conversion of SiO into Si02 occurs in a 
region which is too cold and optically thick to contribute to 
the SiO line intensities, and so the uncertainty in the rate 
coefficient for reaction [6] is not significant in the present 
context. 

Re-adsorption of molecules on to grains ("freeze-out") 
in the postshock gas is included, as in Sch97. The effects of 
freeze-out on the abundance of SiO and its spectrum are 
considered below. 

1.5. Line transfer 

The physical and chemical profiles which derive from the 
shock model summarized above are used in a large velocity 
gradient (LVG) calculation of the intensities of the emis- 
sion lines of SiO and of CO. Our implementation of this 
technique is described in Appendix IA1 where we note that 
an expression for the escape probability which differs from 
that of Sch97 has been adopted. 

2. Results 

2.1. Comparison with the calculations o f\Schilke et al. (1997) 



First, we compare the computed shock structure with that 
of Sch97, for a reference C-type shock model, in which 



the preshock density riH = 10 5 cm~ 3 and the magnetic 
field strength B — 200 fiG, and the shock velocity v s — 
30 km s _1 . The most striking difference between the current 
and the previous model is that the width of the shock wave 
decreases by a factor of approximately 4, to 5 x 10 15 cm, 
from 2 x 10 16 cm in the study of Sch97; see Fig. [Qi. This 
difference is attributable to the more accurate treatment 
of the coupling between the neutral fluid and the charged 
grains in the current model and is an indication of the sig- 
nificance of the inertia of the (negatively) charged grains in 
dark clouds, in which the degree of ionization is low. With 
the narrower shock wave is associated a higher maximum 
temperature of the neutral fluid, as there is less time for the 
initial energy flux, p n v 3 /2, associated with the bulk flow, 
to be converted into internal energy of the H2 molecules or 
to be radiated away. Thus, T n w 4000 K here, compared 
with T n w 2000 K in the model of Sch97. 

There are chemical differences between the models also, 
which arc related in part to the changes in the shock struc- 
ture; these differences may be summarized as follows. 

(i) The fraction of the Si in the grain cores which is re- 
leased into the gas phase by sputtering is approxi- 
mately ten times smaller in the current model than 
in the model of Sch97. This reduction is attributable 
partly to the sputtering yields, which have higher 
thresholds and are smaller for olivine (MgFeSiO^ 
than for the amorphous silica (Si02) considered by 
Sch97; but the main reason for the decrease is the 
enhanced coupling between the neutral fluid and the 
charged grains, which reduces the shock width and 
hence the time available to erode the grains. On the 
other hand, the magnitude of the ion-neutral drift 
speed is similar in both calculations. As a consequence 
of the reduction in the shock width, the integrated 
SiO line intensities predicted by the current model 
are smaller, in general, than calculated by Sch97; see 
Section [3l 

(ii) The displacement of the maximum fractional abun- 
dance of SiO, which forms in the gas-phase reac- 
tions ([I]) and ©, from that of Si, which is eroded from 
the grains, is a more significant fraction of the shock 
width in the current model; cf. Fig. [TJd. The initial 
fractional abundance of O2 in the preshock medium 
is lower here, by a factor of approximately 10, than 
in the model of Sch97, delaying the initial formation 
of SiO. The O2 is assumed to be predominantly in 
the form of ice, which is sputtered rapidly from the 
grains in the early stages of development of the shock 
wave, as may be seen from the two orders of mag- 
nitude increase in the fractional abundance of gas- 
phase O2, apparent in Fig. [Tb. The fractional abun- 
dances of O2 and OH decrease subsequently, at high 
kinetic temperatures, owing to their dissociation by 
H in the chemical reactions 02(H, 0)OH and OH(H, 
0)H2. The former reaction, which is endothermic by 
over 8000 K, proves to be less effective in destroying 
O2 over the smaller width of the current shock model 
(see Fig. []J) than was the case in the calculations of 
Sch97. On the other hand, the lower energy threshold 
of 17 K in reaction allows oxidation reactions to 
proceed further into the postshock region, compared 
with Sch97, whose adopted threshold was 111 K. As 
a consequence, conversion of Si into SiO is slower ini- 
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tially but more complete eventually than predicted by 
Sch97. 

SiC>2 is removed more rapidly from the gas phase in 
the current model. The compression is more rapid 
than in the model of Sch97, and so the rate of 
ion of molecules to grains ("freeze-out") is 
hig&er. If the oxidation of SiO in the reaction SiO(OH, 
H)®02 has an activation energy o f several hundred 
l|in ( |Zachariah fc Tsang (1995)} see Section IP]) . 
tncgmaximum fractional abundance of Si02 would be 
reduced still further. 
AlsoJShown in Fig. QJ) are the fractional abundances 
whieiO are obtained neglecting re- adsorption on to the 
grains. The timescale for freeze-out is sufficiently large that 
chemical profiles are modified only in the cold post- 
shock gas, where the flow speed is practically constant and 
the optical depths in the SiO lines are large. Consequently, 
whilst, the effects re-adsorption on the composition of the 
Si^ostghock gas are dramatic, the freeze-out of SiO has no 
---effect on the predicted line intensities. 

Ii chemical equilibrium, the initial fractional abundance 
of Oj: is approximately 10~ 5 , whereas observations with 
the Odin satellite (Pagani et al. 2003, Larsson et al. 2007) 
have-, placed upper limits of n(02)/n(H2) <J 10~ 7 . In view 
of tliese measurements, we have considered two scenar- 
ios, both with an initial gas-phase fractional abundance 
\n(Oi)/nn = 1.0 x 10~ 7 , as a consequence of the freeze- 
out cf oxygen on to grains, but with differing assumptions 
regarding its chemical form in the grain mantles. 



— The molecular oxygen which formed in the gas phase 
was adsorbed on to the grains, where it remained as 
O > ice in the preshock medium, with a fractional abun- 
dance of 1.3 x 10~ 5 , relative to nn- The initial fractional 



Fig. 1. (a) Temperature of the neutral fluid and veloc- 
ity profiles of the neutral and charged fluids, predicted by 
the present C-type shock model. The shock parameters are 
7iH = 10 5 cm -3 and B = 200 /iG in the preshock gas, and 
v a = 30 km s _1 , C — 5 x 10~ 17 s _1 . The fractional gas- 
phase abundances of selected Si-and O-bearing species are 
plotted in panels (b) and (c); cf. Sch97, fig. 2. The broken 
curves in panel (b) are the corresponding results obtained 
when re-adsorption on to grains is neglected. The discon- 
tinuous curves in panel (c) show the effects of assuming the 
initial abundance of O2 ice to be negligible, i.e. the second 
of the two scenarios described in Section 12.11 



al undance of H2O ice is an order of magnitude larger 
than that the fractional abundance of O2 ice. 
A] omic oxygen was adsorbed on to the grains before O2 
was synthesized and subsequently hydrogenated to H 2 
ici) in the grain mantles of the cold preshock medium. 
T: le fractional abundance of H2 O ice increases by only 
25% as a consequence, to 1.3 x 10~ 4 . 

It turns out that the first scenario is practically equiv- 
alent to assuming that the O2 is initially in the gas-phase 
(see Fig. [Tb), as its release from the grain mantles occurs 
early and rapidly (on a timescale of a few years for the 
model in Fig. Q} in the shock wave. On the other hand, the 
second scenario can result in reduced levels of oxidation of 
Si to SiO in the gas-phase (cf . Fig. [1]) , depending on the 
relative importance of reactions [1] and [2] in the oxidation 
process. In what follows, we present results corresponding 
principally to the first scenario, with the second being con- 
sidered mainly in Appendix [Cj 



2.2, A grid of models 

We have computed a grid of shock models with the following 
parameters: 



nn = 10 cm" 
riH = 10 5 cm' 
tt-h = 10 6 cm 



' 3 , v s = 20,25,30,35,40,45,50 km s" 1 ; 
' 3 , v s = 20,25,30,35,40,45 km s" 1 ; 
' 3 , v s = 20,25,27,30,32,34 km s -1 . 
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In all of these models, we characterized the preshock mag- 
netic field strength by 



B = bn H 5 , 

where nu is in cm -3 and B is in /zG; the scaling parameter b 
was taken equal to 1 . (The effect on Si sputtering of varying 
b is discussed in Section[2T3l) The maximum shock speed for 
tih J> 10 4 cm -3 is determined by the collisional dissociation 
of H2, the main coolant, which leads to a thermal runaway 
and a J-discontinuity (Le Bourlot et al. 2002, Flower & 
Pineau des Forets 2003). 

In fact, we computed two grids, one for each of the sce- 
narios concerning the initial distribution of oxygen between 
O2 and H 2 ices, as specified towards the end of the pre- 
vious Section 12.11 We concentrate on the first of these two 
scenarios, but some additional Figures for the second sce- 
nario are given in Appendix [C] as online material. (Our 
results are available in digital tabular format on request to 
the authors.) 

Because of the sharply defined sputtering threshold en- 
ergy of approximately 50 eV, there is negligible sputter- 
ing of Si from the olivine (MgFeSiO.4) for shock speeds of 
20 km s _1 or less. The fractions of the Mg, Si and Fe which 
are released from the olivine into the gas phase are shown 
in Fig. [2] Comparin g Fig. [2] with the corresponding figure 4 
of |May et al. (200"0) whose sputtering yields are used in 
the present calculations, shows that the fractions of Mg, Si 
and Fe which are sputtered from olivine have decreased by 
an order of magnitude. As the same sputtering yields have 
been used in both studies, this change is attributable to the 
reduction in the shock width, resulting from the improved 
treatment of grain-neutral coupling. We note that CO is 
the principal eroding partner (cf. May et al. 2000). 

Fig. [5] shows that the degree of sputtering is, in fact, 
insensitive to the preshock gas density (cf. Caselli et al. 
1997); it depends essentially on the shock speed. Polynomial 
fits of the sputtered fractions of Fe, Si and Mg, as functions 
of the shock speed, are given in Appendix [B] 

In Fig. [31 we display the fractional gas-phase abun- 
dances of Si and SiO, as functions of the relevant flow time. 
Silicon is produced by erosion of the charged grains by col- 
lisions, principally with molecules, at the ion-neutral drift 
speed. Once the drift speed exceeds the sputtering thresh- 
old velocity, the erosion of Si occurs rapidly, as Fig.[3]shows. 
Thus, the flow time which is directly relevant to the release 
of Si into the gas phase is that of the charged fluid, rather 
than that of the neutrals, which is the appropriate measure 
of the total time for formation of SiO. As noted in item (ii) 
of Section [2TT1 there is an additional, chemical delay to the 
conversion of Si into SiO, in reactions (QJ and which 
is apparent in our Fig. [3j owing to the low abundance and 
partial destruction of O2. The magnitude of this delay de- 
pends on the parameters of the model, notably the shock 
speed, v s , and the preshock gas density, tt-h- Conversion is 
almost instantaneous for v s > 30 km s — 1 , rin — 10 6 cm~ 3 , 
when OH is formed abundantly at the start of the shock 
and reaction (3J dominates the oxidation process. 

Fig. 2] shows the variation with shock speed and 
preshock gas density of the fractional abundance of SiO, 
computed through the entire shock wave. It is evident from 
Fig- HI that the duration of the C-type shock wave, as mea- 
sured by the temperature profile, is of the order of 10 4 , 
10 3 and 10 2 yr for preshock gas densities Tin = 10 4 , 10 5 
respectively. The peak SiO abundance is 



reached over similar timescales. It may be seen that the 
highest fractional abundances of SiO are attained for the 
lowest preshock density, njj = 10 4 cm -3 . At higher densi- 
ties, both O2 and OH, which are reactants in ((T|) and @, 
are destroyed by the atomic hydrogen which is produced 
in the shock wave. Thus, the conversion of Si into SiO be- 
comes incomplete at high density, and the gas-phase SiO 
abundance depends on nn even though the Si sputtered 
fraction does not (cf. Fig. [2]). 

2.3. Dependence on the transverse magnetic field strength 

The existence of a magnetic field transverse to the direction 
of propagation is a necessary condition for a C-type shock 
wave to form, and it is instructive to consider the variation 
of the structure of the shock wave with the strength of the 
magnetic field. Energy equipartition arguments, applied to 
the magnetic and thermal energy densities in the preshock 
molecular gas, of particle density n(H 2 ) + n(He) = 0.6nn 
and kinetic temperature T, suggest that B 2 /(8ir) « n^ksT 
and hence that B = bn^ 5 , where b is a scaling parameter 
(cf. Section |2~2")) such that B is in 11G when nn is in cm~ 3 ; 
this is the proportionality adopted in the grid of models 
presented in Section |2~21 In gas of T = 10 K, equipartition 
with the thermal energy implies b = 0.18. However, we note 
that such a low value of b is inconsistent with the existence 
of a steady-state C-type shock wave when nn = 10 5 cm~ 3 
and u s > 10 km s _1 , as the corresponding ion magnetosonic 
speed in the preshock gas (9.7 km s _1 ) is lower than the 
shock speed. 

In Fig. 03 we present results as a function of the trans- 
verse magnetic field strength, for the parameters of the 
model of Sch97: nn = 10 5 cm~ 3 , v s — 30 km s _1 , and a 
magnetic field scaling parameter b in the range 0.5 < b < 5. 
We recall that Sch97 adopted B — 200 uG, corresponding 
to b = 0.63. This value of b is consistent with the analysis of 
Zeeman measurements by Crutcher (1999)| who concluded 
that there was an approximate equipartition of the mag- 
netic and kinetic energy densities in the molecular clouds 
that he had observed. It may be seen from Fig. [5] that in- 
creasing the magnetic field inhibits the release of Si from 
refractory grain cores in the shock wave, owing to the reduc- 
tion in the maximum ion-neutral velocity difference, Av. In 
Section [3. 51 we consider how the magnetic field strength af- 
fects the relative intensities of the rotational transitions of 
SiO. 



3. SiO rotational emission lines 

The observable quantities are the intensities of the rota- 
tional transitions of SiO and the velocity-profiles of these 
emission lines. Having computed the shock structure, we 
evaluate the line intensities and profiles as described in 
Appendix \X\ assuming that the shock is viewed face-on. 



3.1. Physical conditions in the SiO emission region 

Fig. [5] illustrates the variation of physical conditions 
throughout the formation region of the SiO 5-4 rota- 
tional line for our reference model with nn — 10 5 cm~ 3 , 
v s = 30 km s — 1 , and b — 0.63. It may be seen that the line 
is optically thin through most of the hot precursor (where 
the flow speeds of the charged and neutral fluids differ), 
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Fig. 2. The fractions of Mg, Si and Fe, initially in the form 
of olivine (MgFeSiC^), which are released into the gas phase 
by sputtering within a steady-state C-type shock wave. 

due to both the low SiO abundance and the large velocity 
gradient there. Therefore the line intensity is low despite 
a high kinetic temperature. At the rear of the shock wave, 
approaching maximum compression, the synthesis of SiO 
and the steady decrease in velocity gradient eventually raise 
the optical depth in the line, and the 5-4 intensity peaks, 
with the line temperature attaining values close to the lo- 



Fig. 3. The fractional abundances of Si, released into the 
gas phase by the sputtering of olivine (MgFeSiO^, and of 
SiO, which subsequently forms in reactions 4JXJ) and @. In 
the left-hand panels, the independent variable (abscissa) is 
the flow time of the charged fluid: see text, Section [^1 



cal kinetic temperature of the neutral fluid, T n ; that is, the 
line approaches LTE. The intensity then declines rapidly 
as the gas cools; the decline occurs in 500 yr for the model 
shown here. This behaviour is insensitive to the rate of re- 
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adsorption of SiO on to the grains, which occurs over much 
longer timescales. 

In order to illustrate the dependence of conditions in the 
SiO emission region on the shock parameters, we present, 
in Fig. the important physical quantities, evaluated at 
the peak of the SiO 5~4 line, for all models in our grid. We 
have verified that this is equivalent to computing intensity- 
weighted geometric means of the same quantities over the 
region where the 5-4 line intensity is more than 50% of its 
maximum valueQ Hence, it provides a good indication of 
the mean characteristics of the region producing the peak 
of the emission. 

Fig. [7] shows that the SiO emission peak always oc- 
curs in the cool and dense postshock region: T D w 50 K 
(Fig. [Jji,) and a density of 10-40 times the preshock value, 
close to maximum compression (Fig. [7}d). The SiO frac- 
tional abundance, rt(SiO)/rtn, is also approximately equal 
to its maximum value in the shock wave (compare Fig. [7J; 
with Fig. QJ. The trend to lower SiO abundance at higher 
riH, noted in Section |2~21 is clearly visible here. Finally, the 
"LVG parameter", n(SiO)/(dv z /dz), lies typically in the 
range 10 14 — 10 16 cm" 2 km" 1 s, implying that the 5-4 line 
is optically thick at its peak for most models of our grid. In 
Section 13.51 we compare these physical parameters to val- 
ues inferred previously from LVG analyses of observations, 
assuming a uniform slab which fills the beam. 

3.2. Line profiles and peak line temperatures 

In Fig. [51 we compare the intensity profiles of various ro- 
tational lines, as functions of the flow speed of the neutral 
fluid, expressed in the frame of the preshock gas, for our 
reference model: tih = 10 5 cm" 3 , v s = 30 km s" 1 , and 
b = 0.63. The profiles are seen to be narrow (widths of 1— 
2 km s" 1 ), with similar shapes and peaking within 2 km s 
of v s , as expected for compressed material at the rear of the 
shock wave0Note that the transition 2-1 peaks further into 
the cooling flow, because the lower j— levels are repopulated 
from the higher levels as the temperature falls. 

The general shape and centroid velocities are globally 
similar to those found by Sch97 (their fig. 3b where the 
profiles were plotted in the shock frame), although the dif- 
ferences between the various lines are less significant in 
the present calculations. Also, the emission wing from the 
fast precursor, at the start of the shock wave, is weaker in 
the current models, owing to the delay in SiO formation 
(see Section I2.1[) . making our line profiles narrower than 
in Sch97. Note that including local thermal broadening in 
our profile calculations would not significantly change our 
predicted SiO line width of 1-2 km s— 1, because the SiO 
emission peaks at low temperatures, T n <^ 100 K, where the 

Doppler width is v //cT/44m H < 0.1 km s _1 . 

1 Values evaluated at the peak also differ by less than a factor 
2 from the same parameters evaluated at the "median" point 
where the integrated line intensity, TdV(5 — 4), reaches half of 
its total. 

2 The maximum compression of the postshock relative to the 
preshock gas, \/2v s /va, where v B is the shock speed and va is 
the Alfven speed in the preshock gas, occurs when the magnetic 
pressure in the postshock gas is equal to the initial ram pres- 
sure. It follows that the flow speed in the postshock gas cannot 
fall below u m m = wa/v2 = 1.36 km s" 1 in the shock frame, 
where b is the scaling parameter of the magnetic field, defined 
in Section [231 



In the left column of Figure 02 we show the predicted 
peak temperature of the SiO 5-4 line for the grid of models 
considered in Section [2.2[ as well as the variation with j up of 
the peak brightness temperatures of various lines, relative 
to that of 5-4. The relative intensities have the advantage 
of being independent of the beam filling factor, and thus 
they are comparable directly to observations, without prior 
knowledge of the source size. 

The relative peak intensities are within 20% of unity 
for j up < 7 over a broad range of model parameters 
(v s > 30 km s" 1 and hh < 10 6 cm" 3 ), owing to the large 
opacity and near-LTE excitation conditions. For larger val- 
ues of j U p, the relative intensities are more dependent on 
the shock speed and ~ 1 only when the limiting speed is ap- 
proached. The absolute peak brightness temperature in the 
5-4 line is typically 10-50 K for v s > 30 km s _1 , similar to 
the kinetic temperature in the emission region, but drops 
sharply at lower shock speeds, for which the SiO abun- 
dance (and opacity) is small. The broken curves in Figs. (9ji 
and h are the results obtained assuming that the initial 
abundance of O2 ice is negligible, i.e. the second of the two 
scenarios described in Section |2~T1 

3.3. Integrated line intensities 

In the right column of Fig. 02 are presented the integrated 
intensities (denoted TdV) of the rotational emission lines 
of SiO, relative to the 5-4 line, computed for the grid of 
models considered in Section [2~2l There are significant dif- 
ferences between the relative integrated and peak (T pca k) 
line temperatures (right and left columns, respectively, of 
Fig. 03), owing to systematic variations in linewidth with 
j U p, i.e. in the extent of the region where the line is signif- 
icantly excited. The relative integrated intensities of lines 
with j U p > 7 remain the most sensitive to the shock speed. 
As may be seen in Sch97, the maximum value of TdV oc- 
curs at higher values of the rotational quantum number, 
j U p, for higher 7Jh- Although the shock temperature varies 
only weakly with nn (see Fig. [4j), a higher density enhances 
the rates of collisional excitation of the high-j levels, for any 
given value of the shock speed, v s . In Section [4l we explore 
the usefulness of this effect for constraining the preshock 
density, based on a comparison with actual observations. 

The variation of the absolute integrated intensity of the 
5-4 rotational emission line with the shock parameters is 
shown also in Fig. 03i- The bump in TdV of the 5-4 tran- 
sition, for 30 < v s < 40 km s" 1 and 10 4 < n H < 10 5 cm" 3 
is caused by incomplete O2 destruction in the shock wave, 
resulting in more rapid SiO formation and warmer emission 
zones (cf. Fig. [7]). As the broken curves in Fig. [9ji and h 
show, this "bump" is absent in our second scenario, where 
O2 is never abundant in the gas phase. At higher shock 
speeds, the results from the two scenarios become identi- 
cal, as OH dominates the oxidation of Si in both cases. 

3.4. Influence of viewing angle 

Statistically, there is a low probability that a planar shock 
should happen to be viewed face-on. Accordingly, we have 
explored the effects on the SiO rotational line intensities 
of varying the viewing angle, for the case of our reference 
model, using the formula lA.161 derived in Appendix lAl 
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The variations of Xp C ak(5 — 4) and TdV(5 — 4) with view- 
ing angle, 6, are plotted in the bottom panel of Fig. [TO] 
The peak intensity is almost unaffected by the inclination, 
as the line is already optically thick for a face-on view (see 
Equ. IA.16]) . On the other hand, the velocity projection re- 
duces the line width, and the integrated intensity, TdV, 
declines steadily with increasing viewing angle - by up to 
a factor of 2 at 75°. However, such a variation would be 
difficult to deduce from observations, given the typical un- 
certainties in beam filling factors. 

Panels (a) and (b) of Fig. [TO] illustrate the changes in 
the (relative to 5-4) peak and integrated SiO line tempera- 
tures. Significant changes, compared with viewing face-on, 
are seen only for inclinations greater than 60° from the nor- 
mal, and they affect only the optically thin lines j up < 3 
and j U p > 12. The main change in the curves for the relative 
integrated line temperatures (panel (b) of Fig. [TO]) is that 
the maximum occurs at lower j up as 9 increases; this effect 
could be easily confused with a face-on shock of slightly 
lower preshock density (cf. Fig. [5]). The relative peak tem- 
perature (panel (a)) is even more strongly modified, with 
an upward turn of the curve at j up < 4. The latter charac- 
teristic appears to be the only unambiguous signature of a 
viewing angle > 60° in the context of our one-dimensional 
models. 



3.5. Influence of the transverse magnetic field strength 

In Section 12.31 we have seen that the efficiency of sputter- 
ing Si from grains decreases monotonically with increasing 
magnetic field strength. The effect of varying the scaling 
parameter, 6, on the emergent SiO line intensities is shown 
in Fig. 1 1 1 1 for our reference model. 

In panels (a) and (b) of Fig.[Tl]arc plotted the predicted 
line profiles, peak temperatures, and integrated intensities 
of the SiO 5-4 line, for various values of 6. It may be seen 
that the maximum intensity is reached for intermediate val- 
ues of 0.63 < 6 < 1. At smaller 6, SiO is less abundant: the 
shock wave is narrower and hotter, and so O2 is more read- 
ily destroyed by H, resulting in incomplete oxidation of Si 
into SiO. At larger 6, the SiO emission decreases owing to 
less efficient sputtering of Si (see Fig. [5p) at the lower ion- 
neutral drift speeds. In particular, the predicted intensity 
drops from T p0 ak = 10 K at b — 2 to practically zero at 
6 = 3. 

Panels (c) and (d) of Fig. [TT1 show the relative peak and 
integrated line temperatures, as functions of j up , for vari- 
ous values of b < 2 (curves for b > 2 are not shown as they 
lead to negligible SiO emission). It may be seen that the 
values 0.63 <; b <; 1, which give rise to the strongest 5-4 
emission, also yield the highest relative intensities of lines 
from j up > 7. For example, the intensity of the 11-10 line, 
relative to 5-4, is 3 to 4 times larger than when b = 0.5 or 
b = 1.5. Comparison with Fig. [9] suggests that this depen- 
dence on b may be difficult to distinguish observationally 
from variations in shock speed. A less ambiguous indica- 
tion of the value of b might be obtained from the width of 
the cooling zone, which increases as b 2 , from 10 15 cm for 
b = 0.4 to 2 x 10 16 cm for 6 = 2 and for the parameters of 
our reference model (see Fig. [5^). 



4. Comparisons with observations 

4.1. Rotational line profiles 

As Sch97 first noted, the generic SiO line profile predicted 
by steady planar C-type shock waves, with a peak at high 
velocity (in the postshock gas) and a tail at lower velocity 
(in the accelerating precursor), is reminiscent of the SiO line 
profiles in the L1448 molecular jet (Bachiller et al. 1991). 
Similarly, we note that the reversed shape of SiO profiles in 
the L1157 bowshocks, with a peak at low velocity and a high 
velocity tail (Zhang et al. 1995), could arise if the postshock 
gas is stationary in the cloud frame, i.e. if one observes 
the reverse shock, in which the jet is being decelerated. 
However, in either case, the SiO profiles predicted by our 
models remain narrower than those observed, with widths 
of 0.5-2 km s _1 , as compared to the observed widths of 
5-20 km s" 1 in 3"-10" beams. 

Broader line profiles from steady C-type shocks could 
arise if Si was sputtered not only from grain cores, as as- 
sumed here, but also from SiO-containing grain mantles, 
with lower binding energy. Then, the SiO abundance would 
be much enhanced at intermediate velocities, in the precur- 
sor (see Sch97). However, owing to the steep temperature 
decline across the shock wave, this situation results in large 
variations of the line widths and velocity centroids with the 
emitting rotational level, j up (cf. fig. 5 of Sch97). In fact, 
the observed profiles are very similar from line to line, with 
the (8-7)/ (2-1) intensity ratio showing only modest varia- 
tions with velocity (see, for example, fig. 9 of Nisini et al. 
2007) . These observations suggest that the broad SiO lines 
are not attributable entirely to intrinsic velocity gradients 
through a single, planar C-type shock wave. There may be 
several shock-cooling zones inside the beam, each with a 
narrow intrinsic profile, which appear spread out in radial 
velocity owing to a range of inclination angles or propaga- 
tion speeds, in the observer's frame. This conclusion is sup- 
ported by interferometric observations of LI 157 and L1448 
(Guilloteau et al. 1992; Gueth et al. 1998; Benedettini et 
al. 2007) , which reveal systematic velocity gradients across 
the SiO emitting knots (reminiscent, in some cases, of a 
bowshock geometry; Dutrey et al. 1997) down to 2"-3" res- 
olution; at the distances of L1157 and L1448, the angular 
dimensions of the SiO emitting regions in Fig. [8] for ex- 
ample, are a few tenths of an arcsec. Such complex two- 
dimensional modelling lies outside the scope of the present 
paper. However, we argue in Section 14.31 that we may still 
perform a meaningful comparison of our predicted SiO line 
intensities with observations of knots in outflows, without 
reproducing in detail the line profiles, provided that the 
shock conditions do not vary too much across the beam. 

4.2. Narrow SiO lines near ambient velocity 

In addition to the typically broad SiO line profiles men- 
tioned above, some outflow regions such as NGC1333 and 
L1448 exhibit extremely narrow SiO emission lines, with 
Av w 0.5 km s , near rest velocity (Lefloch et al. 1998; 
Codella et al. 1999; Jimenez-Serra et al. 2004, 2005). The 
corresponding SiO abundance of 10~ n -10~ 10 is two to 
three orders of magnitude smaller than in the broad SiO 
components (Codella et al. 1999; Jimenez-Serra et al. 2005). 
Jimenez-Serra et al. proposed that this feature in L1448 
traces a magnetic precursor, where neutral gas is just be- 
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ginning to accelerate and grain species are starting to be 
released into the gas phase. However, as we now explain, 
detailed multifluid shock models do not support this inter- 
pretation. 

Our SiO line profile calculations show that emission 
from a magnetic precursor does not give rise to a narrow 
feature near the speed of the preshock gas. The line in- 
tensity increases as the neutral fluid is accelerated, heated, 
and enriched in SiO - by orders of magnitude by the time 
that grain-sputtering is complete. Indeed, this deduction 
could have been made already, on the basis of figs. 3 and 
5 of Sch97, which cover the entire velocity range relevant 
to predicting the SiO line profiles. Truncation of the pre- 
cursor when the neutral fluid has been accelerated to only 
v n = 0.5 km s _1 would imply a very finely-tuned shock 
age, a circumstance which appears to us to be improb- 
able. Furthermore, an ion-neutral drift speed of at least 
5 km s _1 is needed to start releasing species from grain 
mantles, where binding energies are a few tenths of an 
eV (Flower and Pineau des Forets 1994), and of at least 
20 km s _1 to start sputtering grain cores (May et al. 2000). 
However, the H 13 CO + line does not show evidence of this 
predicted acceleration: its emission peak is shifted by only 
+0.5 km s _1 from the velocity of the ambient gas, like the 
narrow SiO feature (Jimenez-Serra et al. 2004). 

We believe that a more likely explanation of the nar- 
row SiO feature in L1448 is that it traces Si-enriched post- 
shock material that has been decelerated by and mixed with 
the ambient gas, as proposed originally by Lefloch et al. 
(1998) and Codella et al. (1999) in connection with other 
regions. Given a shock speed v s < 30 km s _1 and the high 
ambient density characteristic of Class protostellar en- 
velopes, deceleration could be achieved readily within the 
L1448 flow age of approximately 3500 yr. The low SiO frac- 
tional abundance would then be a consequence of mixing 
with SiO-poor ambient gas. Alternatively, the narrow fea- 
ture might arise in a reverse C-type shock, where outflow 
material at v < 20 km s _1 is brought almost to rest by 
the much denser ambient medium, and the shock speed is 
too low to produce abundant SiO. Both interpretations are 
consistent with NH3 observations of dense gas in the enve- 
lope of the L1448 protostar, with radial velocity and spatial 
extent similar to that of the narrow SiO feature and signs 
of heating near the path of the fast L1448 jet (Curiel et al. 
1999). 

4.3. SiO line intensities 

If our explanation of SiO profile broadening is correct, one 
could in principle recover the parameters of each individual 
emission zone in the beam by analysing the relative in- 
tensity ratios as functions of velocity. Unfortunately, such 
data are currently quite noisy and not yet available for 
a wide range of values of j up . Furthermore, knowledge of 
the beam-filling factor as a function of velocity would be 
necessary to obtain absolute intensities and remove am- 
biguity in the shock parameters; but this would require 
sub-arcsecond angular resolution, which is not yet avail- 
able. Nevertheless, one may still derive some approximate 
beam-averaged shock properties, if all of the shock compo- 
nents have similar excitation conditions, as is suggested by 
single-dish data, which show the line ratios to be insensi- 
tive to velocity, v. In this case, the observed profile will be 
simply a convolution of the individual, narrow shock pro- 



files with the (unknown) filling-factor, <f>{v). The observed 
absolute TdV is simply that for a single shock, multiplied 
by the total beam filling factor of the SiO-emitting region 
in the beam, / — J cj)(v)dv, as inferred from its overall 
size in single-dish maps. The values of TdV for different 
j up , relative to the 5-4 transition, remain unchanged com- 
pared to a single shock, because / cancels out in the ratios, 
thereby enabling direct comparison with our models. In the 
following, we assume that this situation prevails. 

By way of illustration of the applicability of the shock 
models, we show in Fig.[T2]the relative integrated intensities 
of the rotational transitions of SiO observed in the outflow 
sources L1157 and L1448 (Nisini et al. 2007) and predicted 
by the grid of models, whose parameters are specified; in all 
cases, the magnetic field scaling parameter b = 1. As noted 
in Section 13. 5[ this value of b yields the largest relative in- 
tensities of the high-j lines and hence will yield a lower 
limit to the shock speed required to reproduce the observa- 
tions (except at «h = 10 5 cm' 3 , where high-j excitation 
is a non-monotonic function of i> s ; see Fig. [3f). We assume 
also that the shock wave is viewed face-on; if the true incli- 
nation exceeds 60°, this assumption results in the preshock 
density being slightly underestimated (see Section l3T4")) . The 
models shown as the full curves are those which provide the 
best fits to the observations. In order to illustrate how well 
the shock parameters are constrained, we plot as dashed 
curves "near-miss" models that fit most of the data points, 
or fit all points but do not reproduce the absolute intensity 
of the 5-4 line (see below). 

Fig. [12] demonstrates that steady-state C-type shocks 
with Si-sputtering from grain cores can reproduce success- 
fully the relative integrated intensities of SiO lines in these 
molecular outflows. Furthermore, Fig. [9h shows that the 
models can reproduce also the absolute integrated inten- 
sity of SiO 5-4 with the estimated beam filling factors 
/ w 1 in LI 157 and L1448-R4 and / ~ 1/4 in L1448 
Rl and Bl (cf. Nisini et al. 2007). Alternative fits with 
higher density and lower shock speeds (rin = 10 6 cm -3 and 
27 <; v s <; 30 km s _1 ; n H = 10 5 cm -3 and v s — 25 km s _1 ) 
underestimate TdV (5 — 4) and are thus ruled out. The 
shock speed, v s , is constrained to within 15 km s _1 and 
the preshock density to within a factor of 10, with one of 
our grid values of rijj yielding a clear best fit in all cases. 

In Appendix [Cl we present, in Fig. IC4T results equiva- 
lent to those in Fig.[T2J but for the secondary grid of models, 
in which oxygen is initially in the form of H2O ice rather 
than O2 ice. Again, the relative intensities can be well re- 
produced by steady-state C-type shocks. The absolute SiO 
TdV (5 — 4) favour the low nn, high v s cases (cf. the dashed 
curved in Fig. [9fr) . The best-fit shock parameters and the 
inferred physical conditions at the (5-4) line peak are al- 
most unchanged, as the range of shock speeds is such that 
oxidation of Si by OH is dominant. 

It is instructive to compare the physical parameters at 
the SiO peak of our best-fit, steady-state C-type shock 
models to those previously inferred from an LVG analysis, 
assuming a slab of constant density, temperature, and ve- 
locity gradient along the line of sight (Nisini et al. 2007). 
From Table 1, it may be seen that the 'slab LVG' approach 
yields similar values of the density to our shock models but 
overestimates by a factor 5-10 the kinetic temperature and 
underestimates by several orders of magnitude the Sobolev 
LVG opacity parameter. The cool postshock layer emits 
over a narrow velocity range and needs to be more opti- 
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Table 1. Properties of SiO-emission regions deduced from 
face-on C-type shock models or homogeneous-slab LVG 
models. 



cally thick in SiO to produce the same TdV as a hot slab 
with a large velocity gradient along the line of sight. On 
the other hand, similar SiO abundances are deduced using 
both approaches, to within typically a factor of 3. 

We note that the models should be able to simulate 
also the other spectral observations of the sources, notably 
the H 2 line intensities. In a forthcoming publication, we 
shall consider in detail the outflow L1157 and make a more 
comprehensive comparison of its observed spectrum with 
the predictions of shock models. 

5. Concluding remarks 

We have considered the structure of C-type shock waves 
propagating into molecular gas containing amorphous car- 
bon and silicate grains; olivine (MgFeSiO^ was chosen as 
the representative silicate-grain material. We find that the 
degree of sputtering of silicon from the grains is smaller, 
by about an order of magnitude, than was predicted by 
the calculations of Sch97, owing partly to higher sputtering 
thresholds and lower sputtering yields but principally to 
the reduced width of the shock wave in the present calcu- 
lations. The reduction in the shock width is a consequence 
of a more accurate treatment of the coupling between the 
neutral fluid and the charged grains in the current model. 

A grid of C-type shock models has been computed, for 
values of the preshock gas density and the shock speed 
which are believed to span the ranges of these parameters in 
molecular outflows; two scenarios were considered regard- 
ing the initial distribution of oxygen in the gas and solid 
phases. The maximum speed of a C-type shock is limited by 
the inertia of the (charged) grains in the preshock gas and 
by the collisional dissociation of molecular hydrogen within 
the shock wave, which can lead to a J-type discontinuity 
(Le Bourlot et al. 2002, Flower & Pineau des Forets 2003). 
We find that significant sputtering of the grain cores occurs 
only for shock speeds v s > 25 km s _1 and moderate mag- 
netic field strengths close to equipartition with the cloud 
kinetic energy (magnetic field parameter 0.5 < b < 2). 
However, we note that the sputtering threshold energy is 
determined principally by the so-called "displacement en- 



20 «6 fiel 5 

50-loWe ftp4> th§tgin tlgggabsence of silicon-containing grain 
mantles, SiO line emission in steady-state planar C-type 
shock waves arises predominantly from cool postshock gas, 
close to maximum compression, with negligible emission 
from the precursor. Except at the lowest shock speeds, 
the SiO emission is optically thick and close to LTE for 
4 < j U p < 7. The relative line intensities, as functions of 
j U p, together with the absolute 5-4 line intensity, provide 
good diagnostics of the shock parameters, uh and v s . The 
influence of the viewing angle and transverse magnetic field 
strength is found to be relatively minor over the typical 
ranges of their values. 

Our shock models provide good fits to both the rela- 
tive and absolute SiO intensities in the molecular outflows 
L1157 and L1448; the SiO fractional abundance is deduced 
to be in the range 4 x 10~ 8 < n(SiO)/n H < 3 x 10~ 7 . The 
emission regions of the shock wave are much colder, more 
optically thick, and have 10 to 100 times greater SiO col- 
umn density than estimated previously from optically thin 
LVG slab models (Nisini et al. 2007). Our results are in line 
with a recent analysis of interferometric maps of the HH212 
jet (Cabrit et al. 2007), demonstrating that the SiO emis- 
sion is optically thick and close to LTE, with an intrinsic 
peak brightness temperature of approximately 50 K. On the 
other hand, the line profiles predicted by our planar C-type 
shocks are typically 10 times narrower than observed in 
L1157 and L1448, suggesting that the single-dish beam in- 
cludes shocks with various inclinations and speeds, and/or 
mixing layers. Detailed modelling of the line ratios as func- 
tions of velocity, with a careful correction for differing beam 
widths, would be needed to clarify this issue. 

Whilst the grid of models presented here is intended 
to provide a guide to interpreting observations of outflow 
sources, it should be recalled that the dynamical timescales 
which characterize these regions are often too short to en- 
able C-type shock waves to attain their steady-state struc- 
ture (Chieze et al. 1998, Lcsaffre et al. 2004). The presence 
of an embedded J-type discontinuity has then to be consid- 
ered when modelling specific sources. Studies of the outflow 
source in Orion (Le Bourlot et al. 2002), of jets associated 
with low-mass star formation (Giannini et al. 2004, 2006; 
McCoey et al. 2004), and of the supernova remnant IC 443 
(Cesarsky et al. 1999) have shown that the rovibrational 
line spectrum of H2 can be reproduced successfully by hy- 
brid shock waves, but not by pure C- or J-type shocks. The 
influence of the discontinuity on the C-component (mag- 
netic precursor), and hence on the formation of SiO and its 
emission line spectrum, becomes important for ages smaller 
than the flow time to maximum compression. 

The existence of Si-containing grain mantles is another 
circumstance that would significantly modify the intensities 
of SiO lines and their profiles. The fractional abundance 
of SiO in the warm shock precursor would be enhanced, 
compared with the models considered here, in which Si is 
sputtered exclusively from olivine grain cores. The effects of 
non-steady C-type shocks and Si-containing mantles will 
be considered in a forthcoming paper. 
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Appendix A: SiO radiative transfer 

A.l. Photon escape probabilities 

The SiO rotational level populations and excitation temper- 
atures in our planar shock models are calculated by means 
of a large velocity gradient (LVG) method. This approx- 
imate treatment assumes that, owing to the macroscopic 
velocity field and resulting Doppler shifts, emitted line pho- 
tons are either re-absorbed locally or escape to infinity. The 
escape probability of a line photon in a given direction s is 
then (see Surdej (1977) for a pedagogical derivation): 



(3 8 = 



1 - e" T » 



where the 'LVG optical depth' t s is defined as 



he 



4-7T d(v.s)/ds 



5, 



in 



g\n u 



(A.l) 



(A.2) 



with d(v.s)/ds the radial velocity gradient along direction 
s, ni and n u the number density of molecules in the lower 
and upper levels of the transition, respectively, g\ and g n 
the corresponding statistical weights, and B\ u the Einstein 
coefficient for stimulated absorption. The mean intensity of 
the radiation field at the local frequency v of the transition, 
averaged over all angles, may then be expressed in terms of 
local quantities only as 



J v = S„(l - p) + Ij, 



(A.3) 



where f3 = j (3 s (KI/Att is the escape probability averaged 
over all solid angles; I c is the mean intensity of the contin- 
uum radiation field, taken to be a blackbody (Planck) func- 
tion B V (T) at the temperature of the cosmic background, 
Tbg = 2.73 K; and S v = B U (T CX ) is the local source func- 
tion, where the excitation temperature T ex of the transition 
is defined through n u /n\ = g u /gi exp(— hv /k^T cx ). Two ex- 
pressions have been used for the average escape probability 
0: 

- that of Neufeld & Kaufman (1993): 



Ppla 



1 



f + 3rj_ 



(A.4) 



where t± is the LVG opacity in the z-direction, normal to 
the shock front. This approximation to (3 is accurate for 
a plane-parallel flow, where d(v.M)/ds = n 2 (dv z /dz) and 
t s = tl//i 2 (with the usual notation /i = cos# = s.z). 
- an isotropic approximation: 



isotropic 



= P± = 



1 



(A.5) 
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The first expression was ultimately adopted in the present 
study for consistency with our one-dimensional shock ge- 
ometry, whereas the second was used by Sch97. Since j3 
is always smaller in the plane parallel case (owing to the 
reduced velocity gradients at small /j,), photon trapping is 
more efficient and the excitation temperatures are increased 
compared to the isotropic approximation. Thus, adopting 
the expression of Neufeld & Kaufman (1993) leads to larger 
integrated SiO line intensities; this effect is illustrated in 
Figure IA.ll for our reference shock model. It is marginally 
significant for small and large values of the rotational quan- 
tum number but reaches a factor 3 for j up « 12. 

A necessary condition for the local LVG approximation 
to be valid is that, over the characteristic distance L where 
physical conditions vary, the velocity shift Av arising from 
the velocity gradient should be larger than the local line 
thermal width u t h- Then, line photons are re- absorbed only 
within a region of size < L, where physical conditions and 
SiO excitation are uniform. This criterion can be rewritten 
as 



\dv z /dz\ > Vth/L, 



(A.6) 



where 



i'th 



T n is the temperature of the neutral fluid, and m is the 
mass of the molecule. 

In Figure [A~2l we compare the left and right hand sides of 
(IA.6P for our reference shock model, using L w T n /\dT n /dz\ 
as a characteristic distance. The LVG criterion may be seen 
to be verified throughout the cooling flow of the shock wave, 
where the bulk of SiO emission arises. It is not verified in the 
far postshock region, where the computed velocity gradient 
tends to zero; but this region makes a negligible contribu- 
tion to the line flux, owing to the low escape probabilities. 

A. 2. Numerical implementation 

In the diatomic SiO molecule, electric dipole transitions 
take place only between adjacent rotational levels of the 
ground vibrational state (Aj = ±1 with j the rotational 
quantum number), while collisionally induced transitions 
can, in principle, connect any pair of levels (in practice they 
become less probable with increasing level separation). The 
evolution of the population density m of level i may thus 
be written in the following matrix form: 

drii - 

— n i (A ii _i + B ii+ iJ ii+ i + Bii_i Jj-iji) 

+ 7li-iBi-i t iJi-i t i 

+ '^2'^2n coi \(n j C j i-n i C l j) 1 (A. 7) 

coll j^i 

where n co n is the number density of each collisional partner 
(here H2 and He), CV,- is the collisional rate coefficient from 
level i to level j, and A and B denote the Einstein coef- 
ficients for spontaneous and induced radiative transitions. 
Two approaches to solving for the level populations may be 
taken: 



shock wave, as is done currently for the H2 molecule 
(Le Bourlot et al. 2002). 
— Assume local statistical equilibrium, i.e. set the lhs of 
Equs. fO|) to zero, and solve a posteriori the resulting 
algebraic equations by matrix inversion, using the phys- 
ical structure provided by the MHD shock code. Since 
the radiative terms J depend indirectly on the level pop- 
ulations through the escape probabilities and excitation 
temperatures (E_qu. (|A.3[) ), the procedure must be iter- 
ated, updating J with values of [3 and T cx obtained from 
the previous iteration. 

The former approach is preferable, particularly when the 
molecule under consideration is an important shock coolant 
(e.g. H2) - which is not the case of SiO. Furthermore, the 
flow time in the SiO emission zone is sufficiently long that 
the assumption of local statistical equilibrium is justified. 
Accordingly, we adopted the latter approach in the present 
work. It is common to most applications of the LVG method 
to astrophysical problems (e.g. Sch97, Neufeld & Kaufman 
1993) and has the advantage of being less demanding in 
CPU time. For more significant coolants, such as CO, 13 CO, 
and H2O, the rate of cooling was computed in parallel 
with the shock dynamics, using the c ooling functions calcu- 
lated by Neufeld & Kaufman (1993) by means of the LVG 
method. 

Note that Equs. (|A~7| may be put in a different (matrix) 
form by expressing explicitly J as a function of i?„(T ex ) 
and I c . Using the standard relationships between Einstein 



coefficients, and the definition of T rx in terms of 



and 



— Integrate the set of Equs. (|A.7|) in parallel with the dy- 
namical, thermal, and chemical rate equations of the 



Tij+i, all terms involving T cx cancel from the equations, 
leaving radiative terms which are proportional to the p"s 
(see Goldreich & Kwan 1974): 

diii — 

-j— = + (Aj+l,i + Bi+i^ilc) 

— rii [fli-ijAij-i + pi,i+iBi t i+iI c + /3-j-i,i-Bi,i-iJ c ) 

+ ^i-lPi-l.i-Bj-l.j-fc 

+ ^2 X! n con(njCji - riiCij). (A. 8) 

coll j^i 

Although the two formulations are equivalent and con- 
verge to the same solution, we have found that inversion 
of the latter matrix encounters numerical instabilities and 
convergence problems at high optical depths, possibly due 
to round-off errors in the (vanishingly small) (3 terms. In 
the first formulation, the radiative elements of the matrix 
are never zero, even at high opacity, and we obtain much 
better convergence and accuracy. Hence we have adopted 
(IA.7P in the present calculations. The "lambda-iteration" 
is terminated when we reach convergence of the level pop- 
ulations to 1 part in 10 4 . 

Our code was tested thoroughly _ against the 
routine used by Sch97, in the case (3= f3±, and 
against the web-based online version of RADEX 
(www.sron.rug.nl/ vdtak/radex/radex.php), which 
uses yet another expression for the escape probability, 
valid for a turbulent homogeneous sphere (Van der Tak et 
al. 2007). Discrepancies never exceeded a few percent. 

The MHD shock code provides the physical and chem- 
ical profiles (of the temperatures, densities, velocities, and 
abundances) which are required in order to apply the 
LVG technique. For SiO, we used th e rate coefficients for 
collisional de-excitation published by Turner et al. (1992) 
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for which the collision partner is ground state para-El^. 
These data are available for rotational quantum numbers 
< J < 20 and for kinetic temperatures T = 20, 40, 70, 
100, 150, 200, 250, 300 K and are interpolated to interme- 
diate v alues of T. We made use also of the extrapolated 
data of Schoier et al. (2005) which extend to J = 40 and 
T = 2000 K. Subsequent calculation s, using the rate co- 
efficients of Dayou & Balanga (2006) for collisions of SiO 
with para-H2, have shown that the rotational line inten- 
sities are insensitive to these collision rates because the 
lines are formed under conditions which approach LTE. For 
collisions of SiO with He, we used the rate coefficients of 



Dayou & Balanga (2006) for J < 26 and kinetic tempera- 
tures in the range 10 < T < 300 K. At temperatures higher 
than the maximum for which the rate coefficients were cal- 
culated, we assumed that they remain constant. Upwards 
(excitation) rate coefficients were obtained from detailed 
balance. 

The Einstein A-values, and the rotational constant 
B Q = 21711.967 MHz (corresponding to hB /k B = 
1.042 K) for the ground vibrational state of 28 Si 16 were 
taken from the NIST database (www.nist.gov/data/). 

A. 3. Emergent line intensities and radiation temperatures 

Each layer of the shock wave, of thickness \Az\, elementary 
surface AS, and velocity v z , emits the following luminosity 
(in erg s _1 ) over all directions in a given transition j + 1 — > j 
of frequency v = 2B {j + 1): 



F v (v z ) = (3hvn j+ iA ]+ i^ | Az \ AS. 



(A.9) 



In order to compute line profiles, however, we need to eval- 
uate the specific intensity per unit solid angle, frequency in- 
terval, and projected emitting area (in erg cm~ 2 s _1 Hz -1 
sr _1 ). Following Sch97, we assume that the shock is viewed 
face-on, i.e. in the z-direction, normal to the shock front. 
The intensity is 



I±(v z ) 



47T 



Az 



(A.10) 



where \Av z | is the Doppler width of the layer, viewed along 
the z-axis: 



|Ai/J = -\Av z 



\Az\ x \dv z /dz\ 



Consequently the intensity becomes 

I±{Vz) = ^\d^Jte\ nj+lAj+1 > j - 



(A.ll) 



(A.12) 



Noting that /3±— (1 — e~ T± )/r± and using the definition of 
t in Equ. (|A.2|) as well as the relationships between A and 
_B-Einstein coefficients, one obtains the alternative, simpler 
expression 

I±{v z ) = B v (T eK ) (1 - e- T ^) (A.13) 

which corresponds to the well-known radiative transfer re- 
sult for a uniform slab of excitation temperature T ex and 
opacity t± . 

In order to compare with actual observed SiO spectra, 
one needs to add the cosmic background, and take into 
account the 0N-0FF subtraction applied to radio spectra 
(to remove atmospheric noise). The cosmic background at 
the ON position and velocity v z is -B„(Tb g ) exp(-Tj_), due 



to attenuation by the layer, while at the OFF position it is 
simply BvlTbg). The observed intensity is thus 

lf s (v z ) = ON - OFF = [B tf (T«) - B„(T bg )] (1 - e" TJ -) 

= I ± (v z ) x [l-B v (T hs )/B u (Tl$}lA) 

Finally, we convert the observed specific intensity of the 
layer to a line radiation temperature Tr (in kelvin) using 
the definition 

jobs 2 

T R = 777—^ ( A -!5) 



2h B v 2 

which is standard in radioastronomy. Once radiation tem- 
peratures are calculated for each shock layer, and thus each 
v z , the line profile is integrated over v z to yield integrated 
line intensities, TdV, in K km s _1 . Note that the specific 
intensities, radiation temperatures and TdV are valid for a 
shock that entirely fills the observing beam. Otherwise, the 
observed brightness temperature has to be reduced by an 
appropriate "surface filling factor", / ks AS /(beam area). 

A. 4. Line profile at arbitrary inclinations 

In the general case of an arbitrary viewing angle, /i = cos 9, 
the term t± in Equ. (|A.13p is replaced by the LVG opacity 
in the chosen direction, r(p) = T1//1 2 . Therefore 



(l - c- r ^ 2 ) 



(A.16) 



Thus, Tr will be multiplied by l//i 2 in the optically thin 
regime and will be unchanged in the optically thick regime. 
At the same time, the velocity of the layer, v z , is replaced by 
its projection along the photon path, fj/u z , so the line profile 
becomes narrower. A view which is not face-on will leads to 
a larger TdV for optically thin lines, and to a smaller TdV 
for optically thick lines (assuming that the source still fills 
the beam entirely). A quantitative evaluation of this effect, 
for our reference shock model, is presented in Section [ 



Appendix B: Fractions of Fe, Si and Mg sputtered 
from olivine 

In Section [5^1 we presented and discussed the sputtering of 
Fe, Si and Mg from grains composed of olivine (MgFeSiO^. 
In Fig. IB. 11 we show the numerical fits to these results, as 
functions of the shock speed; they are practically indepen- 
dent of the preshock gas density. Also shown in this Figure 
are the limiting speeds of steady-state C-type shock waves 
for preshock densities ny. = 10 5 and 10 6 cm -3 ; the limiting 
speed is lower for the higher density. The limit is associ- 
ated with the thermal runaway which occurs, owing to the 
collisional dissociation of H2 within the shock wave. The 
numerical values of the coefficients of the polynomial fits, 
of the form 

5 



fX) = ^Oi(X)t 



2 = 



are given in Table [BTTl the shock speed, v s , is in km s . 
We have assumed that the magnetic field parameter 6=1. 
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Table B.l. The numerical values of the coefficients, <Zj(X), 
of the polynomial fits to the fractions of X = Fe, Si and Mg 
sputtered from olivine grains. Numbers in parentheses are 
powers of 10. 



X 


a 


Ol 


a-z 


0.3 


04 


Mg 

Si 

Fe 


-0.5322 
-0.4147 
-0.8579 


0.09172 
0.07364 
0.14681 


-0.006035 
-0.004953 
-0.009607 


0.0001879 
0.0001560 
0.0002979 


-2.746(-06) 
-2.279(-06) 
-4.349(-06) 



Appendix C: Initial gas-phase abundance of O2 

We present, in Figs. IC.lHC.4l the results corresponding to 
those in Figs. [41 171 [51 and [T^l respectively, of the main text, 
but derived from our secondary grid of models, in which 
n(02)/nn = 1-0 x 10 -7 initially in the gas phase but the 
excess oxygen is in form of H2O ice, rather than O2 ice 
as in our primary grid. We recall that the initial fractional 
abundance of O2 in chemical equilibrium is 71.(02) /^h ~ 
1CT 5 . 

In these models, O2 is never abundant in the gas phase, 
and Si oxidation occurs only in reaction O with OH. Thus, 
SiO formation is less efficient at intermediate shock speeds 
and not significant at v s — 25 km s _1 . 
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Fig. 4. The fractional abundance of SiO, n(SiO)/riH, com- 
puted for the grid of shock models and plotted as a function 
of distance (left) and neutral flow time (right); n(SiO)/nn 
is negligible for v s <^ 20 km s _1 . In addition, the temper- 
ature of the neutral fluid is plotted. (See also Fig. IC.ll of 
Appendix [Cl) 
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Fig. 5. (a) The ion-neutral velocity difference, Ad = \vi — 
v n \, and (b) the fractions of Mg, Si and Fe eroded from 
olivine (MgFeSiC^) grains, computed as functions of the 
transverse magnetic field strength, B, in the preshock gas; 
n H = 10 5 cm -3 and v s = 30 km s" 1 . Note that B — 6n^ 5 , 
where & is a scaling parameter (cf. Section [^2J) such that B 
is in jiG when nu is in cm -3 . 



Fig. 6. (a) The temperature of the neutral fluid, T n , the 
brightness temperature, T(5 — 4), in the j — 5-4 line, 
and the compression factor, nn/nn (initial); (b) the opti- 
cal depth, T5_4 in the 5-4 transition and the fractional 
abundance of SiO, n(SiO)/nn, as functions of the flow 
time of the neutral fluid, t n . The model parameters are 
n H = 10 5 cm" 3 , v s = 30 km s" 1 , and b = 0.63. 
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Fig. 7. Physical conditions at the position of the peak in 
the SiO 5-4 line intensity [T poa k(5 — 4)] as functions of the 
shock speed, v s , for all models of the grid: (a) neutral tem- 
perature, T n ; (b) total density, riu; (c) the LVG parame- 
ter, n(SiO)/ (dv z /dz), and the fractional abundance of SiO, 
a;(SiO) = n(SiO)/n H . (See also Fig. EH of Appendix 0) 



Fig. 8. The velocity profiles of transitions from rotational 
levels j —> j — 1 of SiO, computed for our reference model, 
viewed face-on. Only those lines detectable from the ground 
are shown. The model parameters are nn = 10 5 cm~ 3 , 
B = 200 /xG, and v s = 30 km s _1 . The flow speed is in 
the reference frame of the preshock gas. The neutral tem- 
perature profile, T n , is shown also, as are indicative values 
of the flow time of the neutral fluid. (See also Fig. IC.3I of 
Appendix ICl) 
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SiO line peak temperature relative to 5-4 SiO integrated intensity rele 





Fig. 9. (a-c) The peak line temperatures, T pea k, of the rota- 
tional emission lines of SiO, relative to the 5-4 line, as func- 
tions of the rotational quantum number of the upper level of 
the transition, j up , for the grid of models in Section l^T^l The 
value of density, «h, of the preshock gas is indicated in each 
panel, (d) The absolute peak brightness temperature of the 
5-4 line, T pea k(5— 4), as a function of shock speed, v s , for all 
three values of the preshock gas density, nn ■ Shock speeds 
in excess of 34 km s _1 are absent when tir- = 10 6 cm' 3 , as 
they give rise to a J-type discontinuity, and the shock wave 
is no longer C-type; see Section [^1 The right-hand panels 
show the corresponding values of the integrated line inten- 
sities, TdV. The values of TdV(5 - 4) observed in L1157 
and L1448 are indicated. In panels (d) and (h), the broken 
curves show the results obtained assumine that the initial 



Fig. 10. Effect of the inclination angle, 8, on (a) the peak, 
and (b) the integrated intensities of the rotational emission 
lines of SiO, relative to the 5-4 line; (c) the integrated and 
peak intensities of the 5-4 line. 
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Fig. 11. (a) The SiO 5-4 line temperature, as a function of 
the flow speed of the neutral fluid, in the reference frame 
of the preshock gas, for the specified values of the magnetic 
field parameter, 6; (b) the peak and integrated intensities of 
the 5-4 line, as functions of b; (c) the integrated and (d) the 
peak intensities of rotational emission lines j up — > j up — 1 of 
SiO, relative to the 5-4 transition, for the specified values of 
b. All calculations for v s = 30 km s _1 and nn = 10 5 cm~ 3 . 



Fig. 12. The relative intensities of the rotational emission 
lines of SiO observed in the outflow sources L1157 and 
L1448 (Nisini et al. 2007: points with error bars) and pre- 
dicted by the C-type shock models (curves) with the param- 
eters (tt-h, v s ) indicated; see text, Section^ The data points 
which are plotted include a correction for differing beam 
sizes, which is significant for low-j lines. Full curves denote 
the best-fitting models of our grid. Broken curves show 
"near-miss" models with a different value of nn, which ei- 
ther yield a worse fit to the data points or do not repro- 
duce the absolute TdV of the 5-4 line. (See also Fig. IC.4I 
of Appendix ICll 
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Fig. A.l. The integrated intensities of the rotational Fig. A. 2. Consistency of the LVG method: velocity gra- 
transitions of SiO. Full curves correspond to the dicnt criterion (|A.6[) . The model parameters are rtn = 



Neufeld & Kaufman (1993) expression (|A.4|) for the escape 10 5 cm~ 3 , B = 200 [iG, and v s = 30 km s _1 
probability in a plane-parallel flow, dashed curves to the 
isotropic case, equation (|A.5[) . The model parameters are 
n H = 10 5 cm" 3 , B = 200 ^G, and v s = 30 km s" 1 . 
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Fig. B.l. Fits of the fractions of Fe, Si and Mg sputtered 
from olivine grains, as functions of the shock speed. The 
points on the curves indicate the limiting speeds of steady- 
state C-type shock waves, for preshock gas densities nu = 
10 5 and 10 6 cm" 3 , with the lower speed corresponding to 
the higher density. The parameters of the fits to the curves 
are given in Table fETTl 



Fig. C.l. As FiglH but assuming that the initial abundance 
of O2 ice is negligible (the second of the two scenarios de- 
scribed in Section I^T|) . 
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Fig. C.2. As Fig[7J but assuming that the initial abundance 
of O2 ice is negligible (the second of the two scenarios de- 
scribed in Section l27Tj) . 



Fig. C.3. As Fig[Sl but assuming that the initial abundance 
of O2 ice is negligible (the second of the two scenarios de- 
scribed in Section mil)- 




Fig. C.4. As Fig Q21 but assuming that the initial abun- 
dance of O2 ice is negligible (the second of the two scenarios 
described in Section E?T|) . 



